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Abstract: It is well known that a D-string ending on a D3, D5 or D7 brane is 
described in terms of a non-commutative fuzzy funnel geometry. In this article, we 
give a numerical study of the fluctuations about this leading geometry. This allows us 
to investigate issues related to the stability and moduli space of these solutions. We 
comment on the comparison to the linearized fluctuations in supergravity. 
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1. Introduction 

Fuzzy spheres are fascinating examples of non-commutative geometry. They have 
obvious application in string theory, where they arise as solutions to matrix brane 
actions[l][2][3][4][5], and may also play a role in a spacetime explanation of the stringy 
exclusion principle [6]. In this article we will be interested in describing small fluctua- 
tions about the so-called "fuzzy funnel" solutions describing D-strings ending on D3[l], 
D5[2] or D7[4] branes. The funnels have the topology R x S k with S k a fuzzy sphere. 
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Although all of the funnels that we consider are time independent, the equations for 
the funnel also admit interesting time dependent solutions[7]. 

An important difference between a point particle and an extended object is the 
fact that an extended object has a set of low energy excitations which arise from small 
vibrations about their equilibrium configuration. These modes have already played an 
important role in the study of black holes in string theory where they account for both 
the microscopic entropy[8] and the low energy Hawking radiation of the black holes [9]. 
There are a number of reasons why a study of the spectrum of small fluctuations of 
the fuzzy funnels is interesting. First, the prescence of tachyonic modes signals an 
instability in the solution. The fuzzy funnel relevant for D-strings ending on D7-branes 
has been given a description from both the point of view of the D-strings and of the D7- 
branes[4]. It turns out that the two descriptions are in good agreement, at leading order 
in N. From the point of view of the D7 brane description, the funnel is described by a 
7+1 dimensional non-Abelian Born-Infeld theory with a specific non-zero field strength 
switched on. It appears [4] that it is possible to lower the energy of the funnel solution, 
by changing this field strength. Thus there are good reasons to question the stability 
of the D1_LD7 funnel. Our analysis shows that the D1_LD7 system is perturbatively 
stable. Second, if there is a family of funnel solutions this will show up as zero modes 
in the small fluctuation spectrum. Third, the geometry of the fuzzy funnel involves a 
fuzzy sphere. For the fuzzy S 4 and S 6 there are new degrees of freedom[10][ll] that 
have no counterpart in the commutative S 4 and S 6 . These degrees of freedom are most 
important at finite N corresponding to situations in which quantum gravity corrections 
can not be neglected. Understanding the role of these new degrees of freedom is an 
important problem. By studying the small fluctuations of the funnel, we will gain some 
insight into these degrees of freedom. 

The fluctuations with which we concern ourselves arise from the motion of the 
branes in spacetime. Consequently we do not consider the fermions or the gauge field 
that resides on the branes worldvolume. To carry out the small fluctuation analysis, 
we need a basis in which to expand our fluctuations. If we were dealing with the ge- 
ometry of the usual commutative sphere, we would simply use the spherical harmonics. 
Since we deal with the fuzzy sphere we will have to use a deformed algebra, where 
the deformation allows for a cut off on the angular momentum. We will see that the 
usual approach to constructing this basis for the fluctuations is frustrated by the non- 
commuting nature of the coordinates. The natural guess - to take traceless symmetric 
combinations of the underlying coordinates - does not lead to a closed algebra. For- 
tunately, Ramgoolam[10] has found an elegant solution to this problem. For the S 2n , 
instead of simply restricting to the traceless symmetric representations, one considers 
all possible representations of SO(2n+ 1), which may conveniently be labelled by Young 
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Diagrams. The cut off on angufar momentum is implemented as a cut off on the number 
of columns of the Young diagrams. Since there are only a finite number of independent 
basis functions in the deformed algebra the problem lends itself to a numerical study. 
In this article we provide this numerical analysis of the small fluctuations for both the 
D1±D5 and the D1±D7 systems. 

In [2] small fluctuations of the D1_LD5 system were studied. The funnel analysis 
was performed where the funnel is narrow and hence it should be reliably described in 
terms of iV D-strings. These results did not match results obtained for a test-string in 
the background of n D5 branes. Since both the D-string and the test string descriptions 
should be valid for large a this disagreement is rather disturbing. One of our goals in 
this article is to resolve this puzzle. Using a Newtonian approximation we are able 
to produce a model which uses the full funnel geometry as opposed to [2] where the 
funnel was effectively replaced by n D5 branes. We then find that the dynamics of the 
fluctuations depends sensitively on the size of the fluctuation. Further, we show that 
it is possible to reproduce the funnel results using our Newtonian toy model. 

In the next section we briefly review the fuzzy funnel solutions for both the D1_LD5 
and the D1_LD7 systems. In section 3 we discuss in detail how a basis for the fluctu- 
ations can be constructed, allowing for a numerical solution of the small fluctuations 
problem. In section 4, we describe our formulation of the small fluctuations problem 
and in section 5 we give the results of this analysis. The analysis of sections 4 and 5 
is only valid for small field values. In section 6 we compare our results to a gravity 
analysis. In section 7 we present our conclusions and indicate some directions in which 
this work can be extended. 



2. Review of the Fuzzy Funnel Solutions 

As mentioned above, the fuzzy funnel solutions are relevant for the description of D- 
strings ending on either D3, D5 or D7 branes. The funnels have the topology R x S k 
with S k a fuzzy sphere with k = 2 for the D1_LD3 system, k — 4 for the D1_LD5 system 
and k — 6 for the D1_LD7 system. We do not concentrate much on the D1_LD3 system in 
this article. The fuzzy funnels can be described using either a worldvolume description 
of the strings, or a worldvolume description of the branes. The two descriptions are in 
perfect agreement. We will focus on the description in terms of the D-strings. 

The low energy effective action for N D-strings is given by the non-Abelian Born- 
Infeld action [12], [13] 
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= -7\ J d 2 aSTrV- det M , (2.1) 
where 

Q« = <5 ij + 2A[$ i ,$ j ], A = 2tt/ s 2 . 

The symmetrized trace prescription [14] (indicated by STr in the above action) in- 
structs us to symmetrize over all permutations of d a <& % and within the trace 
over the gauge group indices, after expanding the square root. We are using static 
gauge so that the worldsheet coordinates are identified with spacetime coordinates as 
t = x° and a = x 9 . The transverse coordinates are now the non-Abelian scalars $\ 
% = 1, 8. These scalars are NxN matrices transforming in the adjoint representation 
of the U(N) gauge symmetry of the worldsheet theory of the Dls. 
The lowest order equation of motion is 

+ = [&,[&,&]]. 

This equation of motion is valid for small $ J , where the funnel is accurately described 
by a collection of Dls. We are interested in studying the fluctuations <5<3>* about the 
leading configuration $ 

= $ Q + 8&. 

The linearized equation for the fluctuation following from this lowest order equation of 
motion is 

(-d 2 + %)6& = [6&, %}} + [6&, %}} + [^o, 

It is this last equation that we analyze in detail in sections 4 and 5. 

In the remainder of this section, we give the leading configurations <E>g. We start 
by briefly reviewing the construction of fuzzy spheres, and then discuss the funnels 
themselves. In the next section, we will provide a basis[10] in terms of which the 
fluctuations <5$ l can be expanded. 

2.1 Fuzzy Spheres 

To construct the fuzzy S 2m sphere, we need to solve the equation 



S = -Ti / d 2 a STr \ -det 



Vab 

-Xd h ¥ 



-A- 



2m+l 

X i X i = cl, 

i=l 



(2.2) 



with X 1 a matrix, 1 the identity matrix and c a constant. There is a simple construction 
of the matrices X 1 , in terms of the Clifford algebra 



The n-fold tensor product of V (the space on which the P matrices act) is written as 
V® n . The X % are now obtained by taking 

x i = (p ® i ® ... ® i + 1 ® r <g> ... ® 1 + ... + 1 <g> i ® ... <g> P) st . 

The subscript st is to indicate that the above X % are to be restricted to the completely 
symmetric tensor product space for 2m = 2,4, and symmetric and traceless tensor 
product space for 2m = 6. Clearly, the X 1 act in V® n . To prove that the above X 1 do 
indeed provide coordinates for the fuzzy sphere, one shows that Yli^ 1 X l X % commutes 
with the generators of SO(2m+l); then by Schur's lemma, it follows that Yll=i~ 1 X l X % 
is proportional to the identity. 

For properties of the matrices X 1 in the case of the fuzzy S 4 , see [3] and in the case 
of the fuzzy S 6 , see [15] 

2.2 D1/D5 System 

We use iV to denote the number of D-strings in the funnel and n to denote the number 
of D5-branes in the funnel. Five of the eight scalars are excited 



{r,T j } = 2S ij , 



i,j = 1,2,..., 



2m + 1. 



% = R(a)X\ 



i = 



1,2,3,4,5, 



where the function R(cr) is, for small R, 



R(a) 



1 



2V2a 



In this case, the matrices X 1 have dimension 



N 



1 



(n + l)(n + 2)(n + 3), 



6 



and c = n{n + 4). 
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2.3 D1/D7 System 

We use N to denote the number of D-strings in the funnel and n to denote the number 
of D7-branes in the funnel. Seven of the eight scalars are excited 




2 = 1,2,. ..,7, 



R(a) 



1 



2y/3a' 



In this case, the matrices X 1 have dimension 



TV = 



(n + l)(n + 2)(n + 3) 2 (n + 4)(n + 5), 



and c = n(n + 6). 

3. Parametrizing the Fluctuations 

To study the geometry of S 2 embedded in R 3 , we can study the algebra of smooth 
functions over S 2 . This algebra can be taken as a commutative polynomial algebra, i.e 
the polynomial ring R[x l ,x 2 ) x' i } together with the constraint 



In other words we can study R[x l , x 2 , x 3 ]/I where / is the ideal generated by the 
polynomial (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 — c. 

It is then natural to study the geometry of the fuzzy S 2 by studying the algebra of 
functions over it. For some fixed positive integer n we replace x 1 , x 2 , x 3 by three n x n 
matrices X 1 , X 2 , X 3 , which generate the SU(2) Lie algebra 



where C^ k = e^ k and where K depends on the integer n. We still need to impose the 
analog of the constraint (3.1). It turns out that there is a one-to-one correspondence 
between the constraints in the case of the commutative S 2 and the constraints in the 
case of the fuzzy S 2 . This one-to-one correspondence between the constraints can be 
traced back to the fact that the X l 's form a closed Lie algebra among themselves. 
Thus, the basis elements for the non-commutative algebra over the fuzzy S 2 , can be 
obtained by replacing x l — > X 1 in the basis elements of the algebra for the commutative 
S 2 . Note however that we have an upper bound on the number of basis elements in 



(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 



(3.1) 



[X\X j ] = iKC ijk X\ 
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the algebra over the fuzzy S 2 , depending on the dimension of the representation of the 
SU(2) Lie algebra. 

For the case of S 4 and fuzzy S 4 or S 6 and fuzzy S e things are not so simple. We 
discuss in detail the case of the fuzzy S 4 . For the geometry of the commutative 5 4 , the 
algebra of functions is given by R[ 5 ] / I' where I' is the ideal generated by 

the polynomial (re 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 — c. For the fuzzy S 4 , we replace 
x*'s by the X"s given above. The X"s do not form a closed Lie algebra. They can be 
extended to the 5*0(5, 1) Lie algebra by including the new generators 

X H = ^x\X j }. 

The need to include these new elements implies that there are some constraints that 
do not follow from the two sided ideal generated by 

(X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 + (X 5 ) 2 - cl. 
To illustrate this point, note that the relation [2] 

X'i.X 1 - X j X ji = 0, 

follows from 

\ X* ((X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 + (X 5 ) 2 - cl) 
- \ ((X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 + (X 5 ) 2 - cl) X 1 

so that it can be written as an element of the two sided ideal; in contrast, the relation 

X ij X j - 4X* = 0, 

only follows after using the 50(5, 1) Lie algebra. This spoils any possible connection 
between the basis functions in the commutative and non- commutative cases. Thus, 
finding a suitable basis in which to expand the fluctuations of the fuzzy S 4 and fuzzy 
S 6 can not simply be obtained by replacing x % — > X % in the basis for the commutative 
S 4 and 5 6 . 

Fortunately, Ramgoolam[10] has given an elegant solution to this problem. In the 
remainder of this section, we describe the basis provided by[10] for both the fuzzy S 4 
(which is relevant for the D1_LD5 system) and for the fuzzy S 6 (which is relevant for 
the D1_LD7 system). In particular, we will describe our algorithm for generating the 
basis that was used in the numerical computation of the spectrum. 
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3.1 D1/D5 System 

In this section we will describe an 5*0(5) covariant basis for the N x N matrices 5<& l . 
Elements of this basis are labelled by irreducible representations of 50(5) which are 
themselves associated with a unique Young diagram. The Young diagrams labelling 
the representations have at most two rows. To translate a specific Young diagram into 
a basis element, we associate a set of indices to the Young diagram. For example, 
consider the representation associated with the Young diagram with two boxes in the 
first row and one box in the second row. A suitable set of indices would be 



k 


I 


m 





Each index is associated with an X 1 ; for the above example, the basis element is built 
from X k , X 1 and X m . The shape of the Young diagram gives a rule for how the X*s 
should be combined: First we symmetrize over indices in the same row 

x k x l x m -> {X k X l + X l X k )X rn = x k x l x m + x l x k x m , 

and then antisymmetrize over indices in the same column 

j£ rn -\- X^ X k X rn ^ X k X^ X™ X m X^ X k -\- x^ x k x rn x^ x m x k 

Numerically, it is straightforward to implement this symmetrization and then the an- 
tisymmetrization. There is an additional complication because we deal with represen- 
tations of the orthogonal group: the above state can be written as 

X k X^ X m X rn X^ X k -\- x^ x k x rn x^ x rn x k c x a x^ x c 

To get a state which belongs to the advertised irreducible representation, we need to 
restrict ourselves not only to tensors C a bc with the correct symmetrization and anti- 
symmetrization properties, but also, tensors that are traceless. Numerically we im- 
plemented this traceless constraint by subtracting out any projections our basis func- 
tion had on irreducible representations corresponding to Young diagrams containing 
a smaller number of boxes. The cut off on angular momentum becomes a cut off on 
the number of columns in the Young diagrams. One can show [10] that the number of 
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independent states in the irreducible representations included by this prescription is 
precisely equal to 



N 2 



(n + l) 2 (n + 2) 2 (n + 3) 2 
36 



and hence that we do indeed obtain a complete basis in which to express the 5$\ 

Numerically, we performed an exhaustive search over all possible placements of 
indices in the allowed Young diagrams. Of course, many of the operators constructed 
in this way were identically zero. Further, not all of the operators constructed were 
independent, due to the tracelessness constraint. Only the independent operators were 
kept. Our numerical results show clearly that one does indeed construct a complete 
orthonormal basis using the procedure outlined in [10]. 

3.2 D1/D7 System 

In this case, one sums over irreducible representations of 5*0(7) which can be labelled 
by Young diagrams with three rows. Again the cut off is implemented on the columns 
of the Young diagram. In this case, one expects a total of[10] 



basis functions. Again, our numerical results are in perfect agreement with this expec- 
tation. 

4. Problem Formulation 

We are now ready to compute the spectrum of the small fluctuations. As discussed in 
section 2, we consider the lowest order equation of motion. Thus our results will only 
be valid for small $q, corresponding to the large a region of both funnels. In this limit 
both funnels look like a collection of D-strings. 

The problem of computing the small fluctuations essentially reduces to solving the 
eigenvalue problem for the operator L defined by (i = 1, 8) 



Lv^(i) = [v^(j), [0W]] + [«^, [v^(j)M + W, [v {C) (i)\] = \(C)v(°\i). 



The fluctuation is labelled by C. The index % labels the components of a particular 
fluctuation. Thus v^ c \i), fluctuation C, is to be viewed as a vector with 8 components; 
each component is an iVx N matrix. Denote the basis elements described in the previous 
section by (3 A , with A = 1, N 2 . Further, introduce a basis for the fluctuations, given 



N 2 



(n + l) 2 (n + 2) 2 (n + 3) 4 (n + 4) 2 (n + 5) 2 
129600 
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by |C) = \{A, i}). All of the components of the vector are the N x N zero 

matrix, except the i th component, which is f3 A . Thanks to the fact that L has a linear 
action on v^ c \ the eigenvectors of L can be expressed as a linear combination of the 
\{A,i}). Writing 



\{Ai}) = 



p A 5 l8 j 



we define the inner product 



({AOK^'}) = S U 'Tr(p A ip A '). 

We will normalize our states so that ({B,j}\{A,i}) = 5 AB 5 lj . Using the completeness 
relation 

£|{AO>({AOI = i, 

A,i 

where 1 is the 8iV 2 x 8iV 2 identity matrix, we can write 

L\{A,i})=Y,({B,W\{A,i})\{B,j}), 

B,j 

where {{B, j}\L\{A, i}) is an 8iV 2 x 8iV 2 matrix. In terms of the eigenvectors of this 
matrix 



J2({B,j}\L\{A,i})u^(A,i) = X(C)u^(B,j), 

A,i 

the eigenvectors of L can be written as 



,(0 



J2^ C) (A,i)\{A,i}). 



Our approach entails computing the matrix (B, j\L\A, i) and diagonalizing it numeri- 
cally to obtain A(C). After making the ansatz 

6& = f( c \a,T)v( c \i), 
the lowest order equation of motion then reduces to 



{-d 2 r + dj-X(C)R\a))f^(a,r) = 0, 
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which can be dealt with analytically. 

Zero modes (modes with A(C) = 0) correspond to massless fluctuations on the 
worldvolume theory of the D-string. It is natural to interpret these modes as Goldstone 
bosons associated with the spontaneous breaking of translation, rotation and boost 
invariance of the worldvolume theory of the brane. In the large a region where we can 
use (for the D1_LD3 system this is true for all a) 

R{<j) oc — , 
a 

we find that /(cr, r) satsifies 

(-a? + #-£) / (c W) = o. 

Modes with A(C) < correspond to modes with p < 0. Making the ansatz f^ c \a, t) = 
e lUT g(<j), we find the following eigenvalue problem 

w 2 9(°) = {-dl - ^) 9 {a) = Hg(a). 

The eigenvalue problem for this quantum mechanical Hamiltonian has been studied 
previously; recent results include [16], [17]. We will review an argument given in [16] 
which provides a simple illustration of the instability of the D1_LD3 funnel for A(C) < 0. 
Consider 

Set ip\(cr) = \f\il)(\a) so that 
This demonstrates that if [16] 



then H is unbounded from below indicating an infinite number of negative energy 
states. It is possible to argue further that[16] this will occur when \p\ > |. 

The existence of negative energy states (i.e. the existence of bound states for 
the above Hamiltonian) implies that we have modes with u 2 < and which are thus 
tachyonic. For a large enough negative eigenvalue we see that the D1_LD3 system has 



fW(*)\*da-\p\f^da 
J\i>(a)\2da 



2 fW(a)\ 2 da-\p\f^da 
jmaWda 



SWWda 
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an infinite number of tachyonic modes. Thus, instabilities of the D1_LD3 funnel show 
up as negative eigenvalues A(C) 1 . 

These results do not immediately generalize to the D1_LD5 and the D1_LD7 funnels, 
because the bound states have their support at small a which is precisely the region in 
which the funnels differ. For example, the Hamiltonian relevant for the D1_LD5 funnel, 
for small a is 



Although this potential is not as "steep" as for the D1_LD3 system, we still expect bound 
states (signalling tachyons) for large enough \p\. Modes with A(C) > correspond to 
a repulsive potential and hence there can be no bound states. 

Thus, in conclusion modes with A(C) < need to be considered carefully as they 
may indeed signal instabilities. 

5. Results for the D1±D5 System 

The size of the matrix ({B, j}\L\{A,i}) grows rapidly with n. For n = 1 there are a 
total of 16 basis elements and the matrix ({B, j}\L\{A,i}) is a 128x128 dimensional 
matrix; for n = 2 there are a total of 100 basis elements and the matrix ({B, j}\L\{A, i}) 
is an 800x800 dimensional matrix; for n = 3, ({B, j}\L\{A, i}} is a 3200x3200 dimen- 
sional matrix. In this section we will discuss the results obtained for n = 1,2. 



If we had been working on a commutative sphere, then modes with an angular momen- 
tum less than or equal to 1 would be expressed in terms of the scalar or the / = 1 partial 
waves. Here the scalar mode is replaced by the identity matrix and the first partial 
waves by the X 1 . We call fluctuations that can be expressed in terms of the identity 
and X 1 the classical contribution. Since the coordinates do not commute we also have 
contributions that are expressed in terms of [X k , X 1 ]. These have no analog in the case 
of a commutative geometry and hence we refer to these fluctuations as the quantum 
contribution. At this order and at n — 1, to obtain the full set of fluctuations, we have 
checked that there is no need to mix the new quantum degrees of freedom and the de- 
grees of freedom that would be present for a commutative sphere. However, we expect 
this to change once interactions are included and have verified that there is mixing for 

1 The eigenvalues A(C) for the D1_LD3 system are in fact all positive, as expected, since the D1_LD3 
system is supersymmetric and hence obviously stable. 




5.1 71= 1 
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n > 1. Further, looking at the eigenvalues of the fluctuations, we see that the dynamics 
does not discriminate against the new degrees of freedom; they have eigenvalues that 
are comparable to the eigenvalues of the degrees of freedom that would be present for 
a commutative sphere. The spectrum of eigenvalues is 

Eigenvalue Degeneracy Classical Contribution Quantum Contribution 



43 8 35 

8 44 14 30 

16 30 25 5 

24 10 10 

48 1 1 



Since there are no negative eigenvalues, we see that to the order we are working, the 
D1/D5 funnel is stable. 

The modes with eigenvalue equal to zero are giving us information about the moduli 
space (that is, the space of solutions) of the fuzzy sphere. Intuition from the commu- 
tative case suggests that this moduli space is generated by performing translations and 
5*0(1,5) transformations on the fivebrane worldvolume. Looking at the above results, 
we find that there are 8 zero modes corresponding to translations of the X 1 

5¥ oc 1. 

The 5*0(1,5) transformations would include 5 boosts and 10 rotations (symmetries of 
5 4 ). However, we have a remaining 35 zero modes that are all excitations involving 
only the new quantum degrees of freedom [X\ X 7 ]. It has already been argued [10], [11] 
that the non-commutative 5 4 can be viewed as a fuzzy S 2 fiber over an 5 4 base. The 
symmetry of a commutative S 2 is 50(3) and of a commutative 5 4 is 50(5). Assuming 
the symmetries acting on the base and fiber act independently one might expect that 
the 10 rotations are replaced by 3 x 10 = 30 rotations when we consider the non- 
commutative 5 4 . Thus, it seems natural to expect 30 rotations and 5 boosts giving 35 
zero modes for the fuzzy funnel in complete agreement with what we have obtained 
numerically. If this interpretation is correct, we see that the zero mode spectrum 
reflects, in a rather transparent way, the extra dimensions present in the fuzzy sphere 
geometries [10] ,[11]. 

5.2 n = 2 

In this case the generic fluctuation is a linear combination of the new quantum degrees 
of freedom and the degrees of freedom that would be present for a commutative sphere. 
The spectrum of eigenvalues is 
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Eigenvalue Degeneracy 



u 


1 AR 


8 


A A 

AA 


16 


300 


24 


40 


32 


140 


40 


42 


48 


71 


72 


10 


80 


5 



Again, the absence of negative eigenvalues implies that the solution is stable. 

We again find 8 zero modes corresponding to translations of the X 1 . The interpre- 
tation of the remaining 140 zero modes does not seem to be so simple. Assuming our 
interpretation for the zero modes given above is correct, the fact that 140 — 5 = 135 
is not a multiple of 30 suggests that the symmetries acting on the base and fiber no 
longer act independently. Clearly, clarifying the connection between the zero modes 
and symmetries of the non-commuative geometry in general remains an interesting 
open problem. 



Results for the D1_LD7 System 

We will only consider the case n — 1. There are 64 basis elements and the matrix 
({B, j}\L\{A, i}) is a 512 x 512 dimensional matrix. In this case, it is again possible to 
classify the small fluctuations as classical or quantum. Here the quantum contributions 
include both [X\X j ] and 



X ijk = x i X j X k + X i X h X i + X k X i X j - X i X k X j - X j X { X k - X k X j X\ 

This is a direct consequence of the fact that the Young diagrams for SO (7) have three 
rows. The spectrum of eigenvalues is 
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Eigenvalue Degeneracy Classical Contribution Quantum Contribution 



113 8 105 

8 210 210 

16 62 27 35 

24 105 28 77 

48 21 21 

72 1 1 



The absence of negative eigenvalues suggests that this solution is perturbatively stable. 
There are 113 zero modes. We expect 8 of these zero modes correspond to translations. 
The 5*0(1, 7) symmetry of the worldvolume theory of the D7 branes includes 7 boosts. 
The remaining 98 zero modes are presumably associated with symmetries of the fuzzy 
S 6 . 

6. Comparison to Supergravity 

In this section we study the linearized equations for the full Born-Infeld action. For 
the D3_LD1 system, these linearized equations correctly reproduce some of the modifi- 
cations of the supergravity background generated by the D3 branes[18],[19],[l]. For the 
D5_LD1 system a test-string in the background generated by n D5-branes does not agree 
with the linearized equation describing small fluctuations of the D1_LD5 funnel [2]. An 
understanding of this discrepancy is required for a proper interpretation of the fuzzy 
funnel. 

6.1 Linearized Born-Infeld 

Following[2] we consider the fluctuation 

with 1 the N x N identity matrix. Plugging this into the D-string action (in the 
background of the funnel <3>o = RG l / 'y/c\) we obtain 

S=-NT 1 J d 2 a\f (1 + R' 2 )(l - A 2 / 2 ) + A 2 f 2 [l + 4i? 4 /(cA 2 )] 

= ~NT 1 J d*a (h - \\*H^ + l/' 2 ) H(a) = (l + ^) ' . 
Only the quadratic terms are kept as these determine the linearized equation of motion 
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= (Hd 2 T -d 2 a )f(r,a) 

It is natural to ask about the dependence of this linearized equation of motion on 
the specific matrix structure assumed for the fluctuation we consider. Towards this 
end, note that to evaluate the D-string action, we need to compute the terms quadratic 
in 5$ 6 in the determinant 
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It is easy to see that to 0((5$ 6 ) 2 ), the a dependence of the coefficients of both (5$ 6 ) 2 
and (5& 6 ) 2 are unchanged. New terms of the form (<5$ 6 ) 2 and 5$ 6 5$' 6 may be gener- 
ated. This is enough to see that the a dependence of the d 2 5Q 6 and d 2 8& 6 terms in the 
linearized equation of motion are independent of the specific matrix structure assumed 
for the fluctuation we consider. 

6.2 Supergravity Analysis 

The results of the previous section should correspond to the fluctuation equation for 
a test D-string in some supergravity background. Assuming a diagonal background 
metric the test D-string action is [2] 

S = -NT X J d 2 a (Vh - l -\ 2 h^ 2 f 2 + ^A 2 vV 2 ) • 

For n D5 branes 

L 2 

ds 2 = hiayvrj^dafda? + h(a)-^(da 2 + (r 2 dfl 2 3 ), h(a) = 1 + — , 
the resulting linearized equation of motion for / can be written as [2] 

(h 2 d 2 T -dj)f(a,r) = 0, a 2 = a 2 + L 2 . 
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This does not match the result following from expansion of the Born-Infeld action, 
as first observed in [2]. Further, as argued in the previous section, this mismatch is 
independent of the specific matrix structure assumed for the fluctuation of the fuzzy 
funnel description. This disagreement is puzzling as both the D-string and test-string 
descriptions should be valid for large a. The disagreement can be traced back to the fact 
that all three funnel solutions have the same large a behaviour, so that the fluctuations 
from all three funnels look similar to the D3_LD1 system. 

6.3 Newtonian Approximation 

In this section we will build a Newtonian gravity description of the funnel. Our anal- 
ysis provides a toy model description of the full dynamics that illustrates some subtle 
points, explaining the mismatch between the supergravity and linearized Born-Infeld 
descriptions. 

To compare the Newtonian and supergravity descriptions, recall the well known 
fact that for a particle moving with a small velocity, in a region of spacetime well 
approximated by a flat metric, the geodesic equation becomes (we take r to be an 
affine parametrization) 



so that we can compare g 00 to the gravitational potential computed in the Newtonian 
theory. 

A segment of the funnel stretching from a to a + do has a mass (d — 2 for the 
D1_LD3 system and d — 4 for the D1_LD5 system) 



In 9 spatial dimensions, the potential of a point mass drops off as 1/r 7 . Thus, the 
potential set up by the funnel is 



where r is the distance from the observation point (x b s ) to the segment of the funnel 
stretching from a to a + da 




dr 2 



dx a (far 1 dg 00 



dr dr 2 dx l ' 



dm = T d+iy /l + R'(a)R d (a)Q d dcr. 





•^funnel segment (•£ , 0", 0, 0) , 
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x obs = (x°,a o ,6,x d o +2 ,0). 

In what follows we will discuss only the D1_LD5 system. We have checked that exactly 
the same behaviour is exhibited by the D1_LD3 system. Consider the integrand in the 
expression for <fi. We will study this integrand as a function of a for different values of 
x®. The behaviour we will discuss below is qualitatively independent of the details of 
R(a). For this reason we will take R(cr) — - to illustrate our arguments. We will set 
a Q = 2. For x\ = 2, the integrand is shown in figure 1. 



0.09 




I 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 

Figure 1: The integrand for the funnel potential for x & = 2. 



The integral is clearly dominated by the contribution coming from a w 0, that is, the 
potential at this observation point is dominated by the contribution from the D5 brane. 
Decreasing x % = 0.8, the integrand behaviour is shown in figure 2. 
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Figure 2: The integrand for the funnel potential for x % = 0.8. 



We now have two significant contributions - a D5 brane contribution at a ps and 
a contribution from a m 2 corresponding to that portion of the funnel closest to the 
observer, which for large a looks like a D-string contribution. This behaviour is easy 
to understand on physical grounds - as we move closer to the funnel, the local geometry 
has a bigger effect. Decreasing x & = 0.1, the integrand behaviour is shown in figure 3. 




Figure 3: The integrand for the funnel potential for x % Q = 0.1. 
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The integrand is now dominated by a D-string contribution - the contribution from 
the D5 brane at a ~ is negligible. Assuming the funnel reproduces the gravitational 
potential in this third region, we have a convincing explanation of why the fluctuations 
described by the funnel are the same for both the D3_LD1 and the D5_LD1 systems: 
only the geometry of the funnel for large a is relevant, and this is the same for both 
funnels. A comment is in order. In our plots the integrand peaks at a = 2 = a a . 
We chose this value of a Q to produce plots which clearly illustrate our point. In the 
computation of small fluctuations which can be compared to a test string one would 
take a Q » 1 large. 

We have not yet discussed what we mean by a "small fluctuation". We have not 
pursued a suitable definition of what a "small" fluctuation in this Newtonian toy model 
is. Rather, we have been content to show that there is a choice for the scaling of the 
fluctuation that (i) ensures that the funnel predictions are reproduced and (ii) for which 
the fluctuations are small. By reproducing the funnel predictions, we mean we find a 
potential with a dependence. Making the ansatz 

_ j_ 

x d+1 oca n , 

we find that with n = ^ for the D3_LD1 system and with n = ^ for the D5_LD1 
system, the a \ dependence for the potential is obtained. To obtain these values we 
have performed a steepest descent evaluation of the potential assuming that a Q is large. 
The integrand has a peak at a pa a Q . It is pleasing that the fluctuations are indeed 
small in the large a Q region we focus on, showing some degree of self consistency. 

To summarize, our toy model shows that the background corresponding to a stack 
of n D5 branes might not be the appropriate background against which to compare 
the funnel results. Further, we have shown that, for a very specific assumption about 
the a a dependence of the fluctuations, the gravitational results and the funnel results 
are consistent. It would be nice to perform a more complete analysis in supergravity, 
including a careful independent determination of the scaling of x d+1 with a Q , to confirm 
that the two are indeed consistent. 

7. Summary 

In this article we have considered the non-commutative fuzzy funnel geometries that 
provide a description of a D-string ending on a D3, D5 or D7 brane. In particular, 
we have focused on small fluctuations about the leading funnel geometry. A numerical 
treatment of these small fluctuations has been given. The connection between zero 
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modes and symmetries of fuzzy sphere geometries has been discussed. This connection 
is particularly transparent for the D1_LD5 system with n — 1. The comparison of 
the linearized fluctuations of the full Born-Infeld action and the description of these 
fluctuations in supergravity has been considered in the framework of a Newtonian toy 
model. An explanation of previous disagreements between the two descriptions has 
been suggested. 

There are a number of ways in which this work can be extended. The connection 
between zero modes and symmetries of the fuzzy sphere deserves to be studied further. 
This will involve a careful study of the geometry of the fuzzy sphere. The results of [20] 
seem to provide a possible approach for this study. Another interesting direction which 
could be followed would involve improving upon our Newtonian toy model of the funnel. 
The construction of a supergravity model of the funnel as well as a demonstration that 
the D-string description agrees with the supergravity description would be satisfying. 
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